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ABSTRACT 

We have used the transition disk model to fit the simultaneous broad band 
(2 — 500 keV) spectrum of Cygnus X-l from OSSE and Ginga observations. In 
this model, the spectrum is produced by saturated Comptonization within the 
inner region of the accretion disk, where the temperature varies rapidly with 
radius. In an earlier attempt, we demonstrated the viability of this model by 
fitting the data from EXOSAT, XMPC balloon and OSSE observations, though 
these were not made simultaneously. Since the source is known to be variable, 
however, the results of this fit were not conclusive. In addition, since only once 
set of observations was used, the good agreement with the data could have been 
a chance occurrence. Here, we improve considerably upon our earlier analysis 
by considering four sets of simultaneous observations of Cygnus X-l, using an 
empirical model to obtain the disk temperature profile. The vertical structure is 
then obtained using this profile and we show that the analysis is self-consistent. 
We demonstrate conclusively that the transition disk spectrum is a better fit 
to the observations than that predicted by the soft photon Comptonization 
model. In particular, although the transition disk model has only one additional 
parameter, the \ 2 value is reduced and there are no systematic residuals. 
Since the temperature profile is obtained by fitting the data, the unknown 
viscosity mechanism need not be specified. The disk structure can then be 
used to infer the viscosity parameter a, which appears to vary with radius and 
luminosity. This behavior can be understood if a depends intrinsically on the 
local parameters such as density, height and temperature. However, due to 
uncertainties in the radiative transfer, quantitative statements regarding the 
variation of a cannot yet be made. 
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thermal — relativity — stars: Cygnus X-l 



-4- 



1. Introduction 



During the last two decades, the X-ray/7-ray spectrum of Cygnus X-l has been 
extensively studied by, e.g., HEAOl, EXOSAT, Ginga and the Compton GRO. This black 
hole system settles into two distinct long term configurations, namely the so-called soft and 
hard states. In the hard state, the X-ray spectrum is dominated by 2 — 200 keV emission, 
which has been interpreted by some to be the result of the inverse Compton scattering of 
cold external photons by a hot (kT rs 50 keV) low density plasma. This plasma could be 
the innermost region of the accretion disk (Shapiro, Lightman & Eardley 1976 ), or it could 
be a corona overlying a relatively cold disk (Liang & Price 1977; Haardt & Maraschi 1993). 
An alternative model is one in which the emission is a sum of multi-temperature Wien 
peaks from the transition region of the disk (Maraschi & Molendi 1990; Misra & Melia 
1996). All these models predict a similar spectrum which is more or less a power-law with 
an exponential cutoff. 

However, it may now be possible to differentiate between these models by combining 
the data from various high-energy instruments, such as OSSE, Ginga and EXOSAT, and 
thereby studying the broad band nature of this source. The initial steps in this direction 
were taken by Chitnis, Rao & Agrawal (1997) and Gierlinski et al. (1997). They concluded 
that a simple soft-photon Comptonization model does not adequately fit the data. Instead, 
significant systematic residuals were obtained, and it appeared that an additional ad hoc 
spectral component is needed to account for the observations. Gierlinski et al. (1997) have 
also pointed out that the spectrum is photon starved which argues against an overlying 
coronal model in which roughly half of the radiated high-energy photons are reflected back 
into the emitting region by the cold disk. 

Subsequently, Misra et al. (1997) showed that the transition disk model (Misra & 
Melia 1996) appears to fit the broad band data without invoking an additional component. 
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In this picture, the structure of the accretion disk is obtained without assuming a priori 
that the effective optical depth in the vertical direction is much greater than unity. This 
leads to an electron temperature that varies rapidly with radius (Maraschi & Molendi 1990; 
Misra & Melia 1996), and the spectrum differs in a significant way from a simple power-law 
with an exponential cut off. This spectrum fits the data well and is not associated with 
any systematic residuals. These results indicated that the transition disk model may be 
a better fit to the observations than the spectrum expected from a single temperature 
external photon Comptonization. 

However, the earlier analysis by Misra et al. (1997) used non-simultaneous observations 
by EXOSAT, XMPC balloon flights and OSSE. Since the source is known to be variable, 
the results of this analysis were not conclusive. Moreover, since only one set of observations 
was used, it was not clear whether the agreement with the transition model was a chance 
occurrence. An additional complication was the fact that a direct x 2 comparison between 
the transition disk and soft photon Comptonization models could not be made since the 
latter was not fit to the same data. It was shown only that the transition model did not 
produce any systematic residuals whereas the Comptonization model (fit to a different data 
set) did give rise to these. 

In this paper, we use four sets of simultaneous observations of Cygnus X-l made with 
Ginga and OSSE. These data have already been analyzed with an isotropic Comptonization 
model by Gierlinski et al. (1997). Here, we reanalyze the data using the transition disk 
model and compare our results with those of the previous study. Following Misra et al. 
(1997) we fit the data to an empirical model instead of making a direct fit to the more 
complicated exact equations. The best fit temperature profile is then used to determine 
the disk structure. In the final step, the viscosity parameter a is inferred from the disk 
structure. 
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In the next section, we briefly describe the empirical model developed by Misra et al. 
(1997). In §3 the observations and the results of the spectral fit to the transition model 
are compared with those obtained by Gierlinski et al. (1997). In §4 we calculate the disk 
structure from the best fit temperature profile and we discuss the self-consistency of our 
analysis. §5 summarizes the basic results. 

2. The Empirical Model 

To keep the calculation tractable, we use an empirical model developed by Misra 
et al. (1997) to fit the data, rather than attempt a solution of the full equations. We 
emphasis that the difference between the transition disk model and the empirical model is 
in the parametrization. Instead of the physical parameters such as the mass accretion rate, 
the black hole mass and the viscosity parameter a, we use empirical parameters like the 
temperature profile and the absolute normalization described below. 

The gravitational energy released per unit area for each side of the disk is (Shakura & 

Sunyaev 1973) 

_ 3GMM 
F ° ~ 8^^^ J(i?) ' (1) 

where J(-R) = 1 — (Ri/R) 1 ^ 2 and Ri is the radius of the last stable orbit, here taken to be 
Ri = 3R g , with R g = 2GM/c 2 . We assume that the temperature profile of the disk is given 
by 

= kT = 0.1 (-Y J m keV. (2) 

\T / 

where r is the radius of the disk at which 9 = 0.1 J m ~ 0.1 keV. Using this profile and 
assuming a priori (but checked later for self-consistency) that the radiative mechanism is 
saturated Comptonization (i.e., a Wien peak), we obtain the flux at earth (Misra et al. 
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1997): 

f E (E) = KE 2 r°(ey 4 exp(-E/e)r- 2 dr photons/sec/cm 2 /keV , (3) 

Jo 

where E is the photon energy in keV and 

K = 5.59 x 10 27 MD~ 2 , (4) 

where D is the distance to the source. This intrinsic spectrum is a function of K, n, m, and 

To this must be added a reflection component, which has been calculated using the 
appropriate Green's function (White, Lightman, & Zdziarski 1988), under the assumption 
that the reflecting medium is neutral and averaging over incident and observing angles. The 
reflected component depends only on the intrinsic spectrum and the solid angle subtended 
by the reflector, which is conveniently normalized by the factor R re f = Q/2ir. Here R re f = 1 
would correspond to an angle averaged reflection spectrum due to a reflector which suspends 
a solid angle 2-ir with respect to the X-ray source. Since the reflected spectrum depends 
on the incident and observing angles, the real solid angle could be different from the best 
fit Rref which is treated here like a parameter. It should be noted that the shape of the 
reflected spectrum taking into account the angle of observation is similar to the angle 
averaged one, however there could be a difference in normalization (Magdziarz & Zdziarski 
1995). Together with the earlier four parameters, the quantity R re f renders the final 
spectrum (which is used to fit the data) a function of five independent parameters. 

Since the temperature profile of the transition region is obtained in this way, the 
unknown viscosity mechanism does not need to be specified in order for us to obtain the 
disk structure. In principle, the disk equations can be inverted and the viscosity parameter 
(a) can then be obtained as a function of radius. As shown by Misra et al. (1997), however, 
this has only a limited value since the inferred a strongly depends on the details of the 
vertical radiative transfer. Here, we show that this is still the case, but since there are now 



four sets of data, it is possible to infer some useful qualitative trends. We will see that a is 
a function of radius and luminosity, which may be understood in the context of its intrinsic 
dependence on the local parameters such as density, height and temperature. 

3. Observations and Results 

Ginga observed Cygnus X-l on several occasions, one of which was 1991 June 6, during 
the OSSE viewing period 2. Gierlinski et al. (1997) have extracted four Ginga (2 — 30 keV) 
data sets from this period and the corresponding near simultaneous OSSE observations to 
perform their spectral fit. We are using the same set of observations for the present study. 
A systematic error of 1% is included in each of the Ginga channels. 

The OSSE data span the energy range from 50 to 1000 keV, and systematic errors 
are included. These are computed from the uncertainties in the low energy calibration 
and response of the detectors, using in-orbit and prelaunch calibration data. The 
energy-dependent systematic errors are expressed as an uncertainty in the effective area in 
the OSSE response and are added in quadrature to the statistical errors. They are most 
important at the lowest energies, contributing an approximate 3% uncertainty in effective 
area at 50 keV, decreasing to 0.3% at 150 keV and above. The details of these observations 
are given in Table 1 of Gierlinski et al. (1997). 

We fit the four sets of data using the empirical model described in the previous section. 
Following Gierlinski et al. (1997), we consider only data above 3 keV, so that the soft X-ray 
excess does not effect the fitting. The relative normalization between Ginga and OSSE is 
allowed to vary freely within ±15 per cent. In all four cases, the best fit normalization 
falls within this range. The spectrum is absorbed by a line of sight column density using 
the cross section given in Morrison & McCammon (1983). A lower limit of 5 x 10 21 cm~ 2 , 
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corresponding to galactic absorption, is imposed on all the optimizations. This is consistent 
with ASCA measurements by Ebisawa et al. (1996). We model the iron line as a Gaussian, 
whose centroid energy is restricted to lie within the range 6.3 — 6.6 keV and whose line 
width is 0.1 — 0.2 keV, consistent with the ASCA measurements (Ebisawa et al. 1996). 

Table 1 summarizes the results of our spectral fits. The x 2 i n each case is comparable 
to or smaller than those obtained by Gierlinski et al. (1997), who used a single temperature 
isotropic external photon Comptonization model. We show in Figure 1 the unfolded 
spectrum and the residuals for set 1. Unlike the single temperature model, there are no 
strong systematic residuals and hence a second component is not required here. Similar 
results were obtained for the other three data sets. This solidly establishes the results 
indicated earlier by Misra et al. (1997), who used non-simultaneous data. The spectral 
parameters (i.e., n, m and r ) are similar for all four data sets, which suggests that the 
spectral shape does not vary strongly with luminosity, as was also pointed out by Gierlinski 
et al. (1997). 

To check the robustness of these results we have re-fitted the data after fixing n = —3.5, 
which is different from the best fit value n ~ —2.5. The ensuing results are shown in Table 
2. Evidently, the x 2 increased by about ~ 10, but the reduced x 2 is still ~ 1 and thus a 
temperature profile corresponding to this set of parameters cannot be excluded. As we 
note in the next section, these two sets of parameters give rise to different disk structures. 
However, since the A% 2 is large, the latter case should be considered an extreme case. 

4. The transition disk 

Using the best fit empirical temperature profile obtained in the previous section, we 
can now derive the physical structure of the disk, subject to a number of assumptions. 
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Table 1: Spectral parameters for the transition disk model 



Model 


Parameters 


set-1 


set-2 


set-3 


set-4 


component 












abs 


N^ (10 22 cm" 2 ) 


0.5 


0.57 ± 0.22 


1.25 ± .22 


1.1 ± .22 


Gaussian 


E Line (keV) 


6.38 ± 0.1 


6.44 ± 0.1 


6.43 ± 0.11 


6.31 ± 0.21 




a (keV) 


0.2 


0.2 


0.2 


0.2 




^-Line 


9.06 ± 1.6 


10.0 ± 1.7 


5.59 ± 1.0 


4.32 ± 1.3 




EW (eV) 


124 


126 


115 


93.3 


transit 


r 


114.9 ± 15.5 


109.4 ± 10.3 


155.9 ± 26.6 


96.0 ± 17.6 




n 


-2.50 ± 0.16 


-2.62 ± 0.24 


-2.18 ± 0.16 


-2.71± 0.25 




m 


0.87 ± 0.12 


1.03 ± 0.18 


0.62 ± 0.11 


0.98 ± 0.2 




-Tly-e/ 


0.642 ± 0.06 


0.644 ± 0.07 


0.39 ± 0.06 


0.513 ± 0.07 




K 


46.8 ± 0.61 


54.0 ± 0.9 


32.75 ± 0.4 


29.98 ± 0.88 




K 

±i -osse 


0.92 ± .02 


0.89 ± .02 


0.99 ± 0.02 


0.97 ± .04 


X 2 (dof) 




58.96 (75) 


67.33 (75) 


60.59 (75) 


63.52 (75) 



1 The normalization constants for the Gaussian line is in units of 10 3 photons cm 2 s l . 
K osse : Normalization of the OSSE spectrum relative to the Ginga spectrum. 
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Table 2: Spectral parameters for the transition disk model with frozen n = —3.5 



Model 


Parameters 


set-1 


set-2 


set-3 


set-4 


component 












abs 


N^ (10 22 cm" 2 ) 


1.18 ± 0.10 


1.09 ± 0.1 


2.27 ± 0.1 


1.55 ± 0.10 


Gaussian 


E Line (keV) 


6.49 ± 0.14 


6.56 ± 0.15 


6.6 


6.39 ± 0.17 




a (keV) 


0.2 


0.2 


0.2 


0.2 




^-Line 


6.5 ± 1.3 


7.9 ± 1.5 


3.0 ± 0.8 


3.1 ±0.8 




EW (eV) 


88 


97.5 


60 


65 


transit 


r a 


63.43 ± 1.8 


67.13 ± 2.1 


63.35 ± 1.9 


64.2 ± 3.7 




n 


-3.5 


-3.5 


-3.5 


-3.5 




m 


1.54 ± 0.07 


1.69 ± 0.08 


1.51 ± 0.07 


1.59 ± 0.14 




JX r ef 


0.849 ± 0.04 


0.807 ± 0.04 


0.702 ± 0.03 


0.62 ± 0.06 




K 


50.2 ± 0.6 


57.5 ± 0.7 


35.38 ± 0.4 


32.19 ± 0.68 




K 


0.91 ± .02 


0.87 ± .02 


0.96 ± 0.02 


0.93 ± .04 


X 2 (dof) 




71.85 (76) 


72.99 (76) 


83.54 (76) 


67.65 (76) 



x The normalization constants for the Gaussian line is in units of 10 3 photons cm 2 s l . 
K osse : Normalization of the OSSE spectrum relative to the Ginga spectrum. 
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The procedure was developed in detail by Misra et al. (1997), so we provide only a brief 
summary here. 

The assumptions that will be checked for self-consistency are 1) that the plasma cools 
via saturated bremsstrahlung self-Comptonization and 2) that advection is not important. 
The structure equations are those for hydrostatic equilibrium, the radiative transfer and 
radiative cooling. These equations, together with Equation (1), completely determine 
the disk structure with the use of the best fit temperature profile. Note that the torque 
equation generally used in standard disk solutions is not needed here, since the necessary 
information is contained in the temperature profile. Thus the unknown viscosity process 
does not need to be specified in this scheme. The handling of the radiative transfer is the 
most uncertain among these equations and is accounted for with the expression 

F 9 = B^ , (5) 

where F g is the energy flux, P r is the radiation pressure and r is the scattering optical 
depth. The uncertainty is represented by the factor B l7 which is taken to be of order unity. 
In the radiative cooling equation, the bremsstrahlung rate is multiplied by the Compton 
Amplification factor (Svensson 1983). Using the proton electron Coulomb energy transfer 
we take into account the possibility that the proton temperature may be different from 
that of the electrons. The total pressure is taken to be the sum of the radiation and gas 
pressures. 

The normalization parameter K (Eq. 4) may be directly related to the mass accretion 
rate (Misra et al. 1997): 

M = 5.3 x 10 17 (—] I — - — ) grams/sec . (6) 

V50/ \2.5kpcy 

We take the distance to the source to be 2.5 kpc and assume that the black hole mass is ten 
solar masses. 
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In Figure 2, we show the variation of temperature, scattering optical depth and pressure 
with radius for set 1, in which B\ = 0.25. These disk characteristics are qualitatively similar 
to those inferred earlier by Misra et al. (1997). In Figure 3, the effective optical depth 
T e ff is plotted as a function of radius (appearing in the lower right hand corner). In the 
same plot, the Compton y parameter y = (4/cT e /m e c 2 )r 2 is shown as a function of radius. 
The fact that r e // < 1 and y ~ 1 justifies the assumption that the radiative mechanism 
is Comptonized bremsstrahlung. To check the self-consistency of our second assumption 
(regarding the insignificance of advection), we have also calculated the ratio of the pressure 
to the number density times the virial temperature (~ m p c 2 /r): 

*/* = £?. < 7 > 

and have plotted it in Figure 3. Since this ratio is always much greater than one, advection 
is not important and the disk is geometrically thin. Analogous self-consistency checks have 
been carried out for the other three observational sets as well. 

Following Misra et al. (1997), we have next calculated the viscosity parameter a 
as a function of radius from the inverted disk equations. Figure 4 shows this quantity 
corresponding to the first data set for two values of B\ and n. Clearly, a depends rather 
strongly on the actual value of Bi, and hence the radiative transfer, as may be inferred from 
a comparison of the solid and dotted curves in this figure. It also depends in an important 
way on n, as may be seen by inspecting the dashed curve in Figure 4. Since the data do not 
discriminate strongly among the variations of a with radius, it is difficult to make concrete 
statements about the viscosity based on these results. 

Nonetheless, with four sets of data corresponding to different luminosities, it is still 
possible to discern qualitative trends. Figure 5 shows the variation of a with radius for 
the best fit temperature profile in each of the four data sets (solid curves). The viscosity 
parameter a seems to vary inversely with the luminosity, which itself scales with the 
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accretion rate. However, a should in principle depend only on local parameters such as 
height, density and temperature. We have found that the following empirical "fit" for a in 
terms of these parameters, reproduces the variation quite well: 

-0.85 / - \ -2.8 



/ T 
a(E,T) =2.65 ' 



V10 8 K/ yiOgcm- 2 ; ' v ' 

where E is the column density. This empirical fit is shown (dotted curve) in Figure 5. Care 
should be taken when interpreting this result since a depends sensitively on the uncertainty 
in the radiative transfer. Moreover, the fit itself may not be unique. We are showing this 
here primarily to demonstrate that it is possible to obtain a functional form for a that 
depends only on the local parameters. 

5. Summary 

We have produced a fit to four sets of simultaneous broad band data for Cygnus X-l 
within the context of the transition disk model. These data were obtained with observations 
by OSSE (50 — 500 keV) and Ginga (3 — 30 keV). We have established rather strongly 
that the transition disk model fits the data better than a pure soft-photon Comptonization 
model, as suggested earlier by Misra et al. (1997) based on non-simultaneous data. In 
particular the \ 2 value for the transition disk model is in each case smaller than or 
comparable to that of the isotropic Comptonization model fit to the same data (Gierlinski 
et al. 1997). In addition, unlike the soft-photon Comptonization model, there are no 
systematic residuals for the transition disk fit. The transition disk is described by four 
parameters as compared to the isotropic Comptonization model which has three. However, 
for the model presented here no additional ad hoc component, with two extra parameters 
(Gierlinski et al. 1997) are needed. 

We have also determined the disk structure from the best fit temperature profiles, and 
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through self-consistency checks, we showed that our initial assumptions were valid. The 
viscosity parameter a inferred from the calculated disk structure was found to depend on 
the details of the vertical radiative transfer and the spectral fitting. Thus, it is not yet 
possible to make concrete statements regarding it's value or specific tendencies. However, 
since we were able to analyze the data at four different epochs (each corresponding to a 
different luminosity), we derived the qualitative variation of a with radius and luminosity, 
which we accounted for in terms of its dependence on the local disk structure parameters, 
such as density, temperature and height. 

The authors thank A. A. Zdziarski and M. Gierlinski for providing us with the reduced 
data from OSSE and Ginga. RM and VRC thank A. Kembhavi for useful discussions. VRC 
acknowledges IUCAA for its hospitality. This work was partially supported by NASA grant 
NAG 5-3075. 
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Fig. 1. — The deconvolved spectrum of Cygnus X-l corresponding to the first data set, 
obtained from simultaneous observations with Ginga and OSSE. The contributions from 
individual model components, viz., the iron line and reflection, are shown separately. The 
residuals to the model fit are shown in the lower panel as a contribution to the x 2 . 

Fig. 2. — (a) The temperature profile in the disk for the best fit values in data set 1, and 
B\ = 0.25. Solid curve: electron temperature; dotted curve: proton temperature, (b) 
The scattering optical depth (r = HpUxH) as a function of radius corresponding to the 
temperature profile in (a) and B\ = 0.25. (c) The mass density as a function of radius 
corresponding to the temperature profile in (a) and B\ = 0.25. (d) Pressure as a function of 
radius corresponding to the temperature profile in (a) and B\ = 0.25. Solid curve: Radiation 
pressure; Dotted curve: Gas pressure. 

Fig. 3. — The effective optical depth (r e //), the Compton y parameter and the ratio P a /P 
as functions of radius for data set 1. 

Fig. 4. — The viscosity parameter a as a function of radius for data set 1. Solid curve: 
B\ = 0.25; Dashed curve: B\ = 0.5. The dotted curve is for a temperature profile calculated 
with n = —3.5. 

Fig. 5. — The viscosity parameter a (solid curve) as a function of radius for data sets 1, 2, 
3, and 4. Dotted curve: empirical fit given by Equation (8) in the text. 
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